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Abstract.  We  present  the  first  full  FIR  spectrum  of  Centaurus 
A  (NGC  5128)  from  43  -  196.7  pm.  The  data  was  obtained  with 
the  ISO  Long  Wavelength  Spectrometer  (LWS).  We  conclude 
that  the  FIR  emission  in  a  70  "  beam  centred  on  the  nucleus 
is  dominated  by  star  formation  rather  than  AGN  activity.  The 
flux  in  the  far-infrared  lines  is  ~  1%  of  the  total  FIR:  the  [C  II] 
line  flux  is  ~  0.4%  FIR  and  the  [O  I]  line  is  ~  0.2%,  with  the 
remainder  arising  from  [O  ill],  [N II]  and  [N  ill]  lines.  These 
are  typical  values  for  starburst  galaxies. 

The  ratio  of  the  [N  III]  /  [N II]  line  intensities  from  the  HII 
regions  in  the  dust  lane  corresponds  to  an  effective  temperature, 
~  35  500  K,  implying  that  the  tip  of  the  main  sequence  is 
headed  by  08.5  stars  and  that  the  starburst  is  ~  6  xlO6  years 
old.  This  suggests  that  the  galaxy  underwent  either  a  recent 
merger  or  a  merger  which  triggered  a  series  of  bursts.  The  N/O 
abundance  ratio  is  consistent  with  the  range  of  ~  0.2  -  0.3  found 
for  Galactic  HII  regions. 

We  estimate  that  <  5%  of  the  observed  [C  II]  arises  in  the 
cold  neutral  medium  (CNM)  and  that  ~  10%  arises  in  the  warm 
ionized  medium  (WIM).  The  main  contributors  to  the  [C  II] 
emission  are  the  PDRs,  which  are  located  throughout  the  dust 
lane  and  in  regions  beyond  where  the  bulk  of  the  molecular  ma¬ 
terial  lies.  On  scales  of  ~  1  kpc  the  average  physical  properties 
of  the  PDRs  are  modelled  with  a  gas  density,  n  ~  103  cm~3, 
an  incident  far-UV  field,  G  ~  102  times  the  local  Galactic  field, 
and  a  gas  temperature  of  ~  250  K. 

Key  words:  galaxies:  individual:  Centaurus  A  =  NGC  5128  - 
infrared:  galaxies  -  galaxies:  ISM  -  galaxies:  starburst  -  galax¬ 
ies:  active 


1.  Introduction 

Centaurus  A  (NGC  5128)  is  the  nearest  (d  =  3.5  Mpc;  1  "~17  pc, 
Hui  et  al.  1993)  example  of  a  giant  elliptical  galaxy  associated 
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with  a  powerful  radio  source.  The  large-scale  radio  morphology 
consists  of  twin  radio  lobes  separated  by  ~  5  degrees  on  the  sky. 
The  compact  (~  milliarcsecond)  radio  nucleus  is  variable  and 
has  a  strong  jet  extending  ~  4  arcminutes  towards  the  northeast 
lobe.  The  spectacular  optical  appearance  is  that  of  a  giant  ellip¬ 
tical  galaxy  that  appears  enveloped  in  a  nearly  edge  on,  warped 
dust  lane.  There  is  also  a  series  of  faint  optical  shells.  The  stellar 
population  in  the  dominant  elliptical  structure  is  old,  whilst  that 
of  the  twisted  dust  lane  is  young,  sporadically  punctuated  by  HII 
regions,  dust  and  gas  (Graham  1979).  The  overall  structure  of 
Cen  A  resembles  that  of  a  recent  (<  4  x  108  years,  Tubbs  1980) 
merger,  between  a  spiral  and  a  large  elliptical  galaxy.  The  dust 
lane  is  the  source  of  most  (90%)  of  the  far-infrared  luminosity 
(Lfir  ~  3  x  109  Lq)  and  is  thought  to  be  re -radiated  starlight 
from  young  stars  in  the  dusty  disk  (Joy  et  al.  1988). 

In  Sect.  2  we  describe  the  observations  and  data  analysis. 
Sect.  3  looks  at  the  general  FIR  properties  and  proceeds  to  model 
the  HII  regions  and  the  PDRs  in  the  dust  lane.  Sect.  4  summarises 
the  results  and  presents  our  conclusions. 

2.  Observations 

Cen  A  was  observed  with  the  LWS  grating  ( R  =  A/AA  ~  200) 
as  part  of  the  LWS  consortium’s  guaranteed  time  extragalactic 
programme.  A  full  grating  observation  (43  -  196.7  pm)  was 
taken  of  the  nucleus  at  the  centre  of  the  dust  lane  and  a  series  of 
line  observations  were  taken  at  two  positions  in  the  SE  and  NW 
regions  of  the  dust  lane.  A  short  [C  II]  157  pm  line  observation 
was  taken  off-source  at  position  #4  (see  Table  1)  to  estimate 
the  Galactic  emission  near  the  source.  Position  #1  was  intended 
to  provide  a  deeper  integration  coincident  with  position  #2,  but 
was  accidently  offset. 

A  series  of  half-second  integration  ramps  were  taken  at 
each  grating  position  with  four  samples  per  resolution  ele¬ 
ment  (AA  =  0.29  pm  AA43  —  93  pm  and  AA  =  0.6  pm 
AA84  —  196  pm).  The  total  integration  time  per  resolution  el¬ 
ement  and  per  pointing  were:  position  #1  88  s  for  the  [O  III] 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2000 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2000  to  00-00-2000 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

ISO-LWS  spectroscopy  of  Centaurus  A:  extended  star  formation 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory ,4555  Overlook  Avenue, 

SW, Washington, DC, 20375 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 

OF  PAGES 

6 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


886 


S  J.  Unger  et  al.:  ISO-LWS  spectroscopy  of  Centaurus  A:  extended  star  formation 


50 


55 


.2  -43:00 
c 

o 

<0 

□ 

05 


10 


15 


Fig.  1.  Cen  A  digital  sky  survey  image  overlaid  with  the  LWS  beam 
positions 


Table  1.  Observation  log 


Position 

Offset  in  RA 

arc sec 

Offset  in  Dec 

arc sec 

Date/AOT 

#1 

+  29 

-  12 

1996  Aug  23  L02 

#2  Centre 

0 

0 

1997  Aug  11  L01 

#3  (map) 
NW 

-53 

+  27 

1997  Feb  12  L02 

SE 

+  110 

-49 

#4  (off) 

-2 

+  600 

1997  Feb  12L02 

Offsets  w.r.t.  13h  25m  27.6s  -43d  01  '08.6  "J2000 


52  /im  and  34  s  for  the  [N  ill]  57  ^im;  position  #2  (the  centre), 
30  s  for  the  range  43-196  /im;  positions  NW  and  SE  (2  point 
raster  map)  22  s  for  the  the  [O  I]  63  /im,  14  s  for  the  [O  ill] 
88  /tm,  12  s  for  the  [N  II]  122  /im,  28  s  for  the  [O  i]  145  /tm 
and  12  s  for  the  [C  II]  158  /im;  position  #4  12  s  for  the  [C  II] 
158  /tm. 

The  data  were  processed  with  RAL  pipeline  7  and  analysed 
using  the  LIA  and  ISAP  packages.  The  LWS  flux  calibration 
and  relative  spectral  response  function  (RSRF)  were  derived 
from  observations  of  Uranus  (Swinyard  et  al.  1998).  The  full 
grating  spectrum  at  the  centre  enabled  us  to  estimate  the  rela¬ 
tive  flux  uncertainty  between  individual  detectors  arising  from 
uncertainties  in  the  relative  responsivity  and  the  dark-current 
subtraction.  The  offsets  between  the  detectors  (excluding  de¬ 
tector  SW1)  was  <  10%.  The  [Om]  88  /im  line  on  detectors 
SW5  and  LW1  had  a  15%  systematic  uncertainty  and  the  [C  II] 
line  on  detectors  LW3  and  LW4  had  a  10%  systematic  uncer¬ 
tainty.  We  therefore  adopt  a  relative  flux  uncertainty  of  ~  15%. 
Because  we  only  took  spectra  of  individual  lines  at  the  NW  and 
SE  positions  there  is  no  corresponding  overlap  in  wavelength 


Wavelenqtn  (/tm ) 

Fig.  2.  LWS  spectrum  of  the  central  region  of  the  dust  lane  (dashed  line 
is  the  blackbody  fit) 


coverage  at  these  positions.  One  indicator  of  relative  flux  uncer¬ 
tainty  is  a  discrete  step  down  in  flux,  of  ~  25%,  at  ~  125  /tm  at 
the  SE  position.  The  relative  flux  uncertainty  is  assumed  to  be 
<  25%  at  these  positions. 

The  absolute  flux  calibration  w.r.t.  Uranus  for  point  like  ob¬ 
jects  observed  on  axis  is  better  than  15%  (Swinyard  et  al.  1998). 
However,  extended  sources  give  rise  either  to  channel  fringes 
or  to  a  spectrum  that  is  not  a  smooth  function  of  wavelength. 
This  is  still  a  calibration  issue.  For  example,  in  Fig.  2,  detec¬ 
tors  SW5,  LW1,  LW2  have  slopes  that  differ  from  those  of  their 
neighbours  in  the  overlap  region.  This  may  account  for  the  con¬ 
tinuum  shape,  which  is  discussed  in  Sect.  3.1.  The  LWS  beam 
profile  is  known  to  be  asymmetric  and  is  still  under  investi¬ 
gation.  We  therefore  adopt  a  value  for  the  FWHM  of  70  "  at 
all  wavelengths,  believing  that  a  more  sophisticated  treatment 
would  not  significantly  affect  our  conclusions.  We  also  note  that 
there  is  good  cross  calibration  between  the  ISO-LWS  results 
and  the  Far-infrared  Imaging  Fabry-Perot  Interferometer  (FIFI) 
(Madden  et  al.  1995);  the  [C  II]  peak  fluxes  agree  to  within  ~ 
10%. 


3.  Results  and  discussion 

3.1.  General  FIR  properties 

The  far-infrared  continuum  at  each  position  is  modelled  with 
a  single-temperature  blackbody  spectrum  of  the  form  Fa  a  fl 
Bv(T)(l-e~rdust),  where  the  solid  angle,  O,  is  constrained  to 
equal  the  LWS  beam,  Ba(T)  is  the  Planck  function  at  temper¬ 
ature  T  and  Tdust  ex  A-15.  The  result  for  the  central  position 
is  shown  as  the  dashed  curve  in  Fig.  2.  Although  the  observed 
continuum  is  not  a  simple  function  of  wavelength  and  the  single 
temperature  blackbody  is  not  an  especially  good  fit,  particularly 
at  wavelengths  >  100/tm,  a  better  calibration  of  straylight  and 
the  beam  profile  is  required  for  anything  more  sophisticated. 
The  best  FIR  temperature  at  each  position  is  ~  30  K. 

The  luminosities  quoted  here  are  derived  from  the  line  fluxes 
listed  in  Table  2.  At  the  central  position,  the  total  luminosity  in 
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Table  2.  Line  fluxes 


Line 

Arest 

/im 

#1 

Flux 

Centre 

Flux 

NW 

Flux 

SE 

Flux 

[O  hi]  3p2-3pl 

51.815 

0.31 

0.72 

- 

- 

[N  ill]  2p|-2pi 

57.3170 

0.16 

0.24 

- 

- 

[Ol]  3pl-3p2 

63.184 

- 

1.96 

0.90 

0.51 

[O  hi]  3pl-3p0 

88.356 

- 

0.70 

0.6 

0.2 

[Nil]  3p2-3pl 

121.898 

- 

0.15 

0.15 

<0.15 

[Ol]  3p0-3pl 

145.525 

- 

0.11 

0.08 

<0.03 

[C  ii]  2pf-2pl 

157.741 

- 

2.91 

2.43 

0.92 

Flux  xlO"18  W  cm-2 
-  wavelength  range  not  covered 

upper  limits  are  3  x  rms  residuals  from  a  fit  to  the  continuum 
[C  II]  157  pm  flux  at  #4  (off)  is  0.04 


all  of  the  far-infrared  lines  is  2.6  x  107  Lq,  which  is  ~  1%  of 
the  total  FIR  luminosity  (L43_1g7Atm  =  3.2  x  109  L0).  The 
[C  II]  luminosity  is  1.1  x  107  L0  (0.4%  FIR)  and  the  [O  i]  lumi¬ 
nosity  is  7.5  x  106  L0  (0.2%  FIR).  Because  full  spectra  are  not 
available  at  the  NW  and  SE  positions  we  estimate  the  FIR  con¬ 
tinuum  luminosity  by  integrating  under  the  single-temperature 
blackbody  fit  to  the  data.  At  the  NW  position  the  total  FIR  lu¬ 
minosity,  L43_1g7Aim  =  2.2  x  109  L0.  The  [C  n]  luminosity  is 
9.3  xlO6  L0  (0.4%  FIR)  and  the  [Ol]  luminosity  is  3.4  xlO6 
L0  (0.2%  FIR).  At  the  SE  position  the  total  FIR  luminosity, 
L43-ig74im  =  8.7  x  108  L0.  The  [C  Ii]  luminosity  is  3.5  xlO6 
L0  (0.4%  FIR)  and  the  [O  I]  luminosity  is  2.0  x  106  L0  (0.2% 
FIR).  These  are  typical  values  for  starburst  galaxies  (c.f.  Lord 
et  al.  1996). 

3.2.  Ionized  gas  lines 

Photons  of  energy  35.12,  29.60  and  14.53  eV  are  required  to 
form  0++,  N++  and  N+,  respectively,  so  that  the  observed 
[O  ill],  [N  ill]  and  [N II]  emission  must  originate  in  or  around 
HII  regions.  The  [O  ill]  line  ratio  is  a  sensitive  function  of  den¬ 
sity  in  the  range  ~  30  —  104  cm'3.  For  the  central  position  this 
ratio  is  ~  0.9,  corresponding  to  an  electron  density,  ne  ~  100 
cm'3  (Rubin  et  al.  1994).  The  [O  III]  lines  indicate  a  higher 
electron  density,  nc  ~  250  cm'3,  for  the  starburst  nuclei  of 
M82  (Colbert  et  al.  1999)  and  M83  (Stacey  et  al.  1999).  In  con¬ 
trast,  the  [N II]  205  pm  /  122  pm  line  intensity  ratio  for  the 
Galaxy  gives  an  average  electron  density,  of  only  ~  3  cm'3 
(Petuchowski  &  Bennett  1993).  The  thermal  pressure  of  the 
ionized  material  in  the  Cen  A  dust  lane  is  therefore  closer  to 
that  of  starburst  galaxies  than  to  that  of  the  Milky  Way. 

Since  the  [N  ill]  57  pm  and  the  [N  n]  1 22  pm  lines  arise  from 
different  ionization  states  of  the  same  element,  the  line  intensity 
ratio  is  sensitive  to  the  hardness  of  the  interstellar  UV  field  and 
therefore  to  the  spectral  type  of  the  hottest  main  sequence  star. 
For  the  central  position  [N  ill]  /  [N  II]  ~  1.6.  This  is  larger  than 
the  value  of  ~  0.9  for  M83  (Stacey  et  al.  1999)  but  smaller  than 
the  value  of  ~  2.1  for  M82  (Colbert  et  al.  1999).  Assuming 
that  the  region  is  ionization  bounded,  with  an  electron  density, 
ne  ~  100  cm'3  the  [N  ill]  /  [N II]  line  intensity  ratio  for  Cen 


A  corresponds  to  an  abundance  ratio  N++/N+  of  ~  0.3;  this 
corresponds  to  an  effective  temperature,  Tcff  ~  35  500  K  (Rubin 
et  al.  1994).  Applying  the  same  corrections  to  M82  and  M83 
with  ne  ~  250  cm'3  implies  an  effective  temperature,  Teg-  ~ 
34  500  K  for  M83  and  Tcff  ~  35  500  K  for  M82.  If  the  effective 
temperature  in  Cen  A  corresponds  to  the  tip  of  the  main  sequence 
formed  in  a  single  starburst,  we  are  observing  08.5  stars  and 
the  burst  is  ~  6  xlO6  years  old.  If  the  burst  was  triggered  by 
the  spiral-elliptical  galaxy  merger  then  its  occurance  was  very 
recent.  Alternatively,  the  merger  triggered  a  series  of  bursts  of 
star  formation,  of  which  we  are  witnessing  the  most  recent. 

The  N++  and  0++  coexist  in  roughly  the  same  ionization 
zones,  and  the  [O  ill]  52  pm  and  [N  ill]  57  pm  lines  have  roughly 
the  same  critical  density.  As  a  result  the  ratio  of  these  lines  is  an 
indicator,  to  within  ~  50%,  of  the  N++/0++  abundance  ratio, 
which  itself,  is  nearly  equal  to  the  N/O  ratio  in  the  hard  UV  field 
environments  we  are  seeing  here  (Rubin  et  al.  1988).  The  line 
ratio  we  observe  at  the  centre  of  the  dust  lane  is  ~  0.3  -  the  same 
as  found  in  the  nucleus  of  M82  (Colbert  et  al.  1999),  but  much 
smaller  than  that  found  for  the  nucleus  of  M83  (~  0.67  Stacey 
et  al.  1999). 

A  more  precise  determination  of  the  abundance  ratio  re¬ 
quires  the  observed  line  ratio  to  be  divided  by  the  volume  emis- 
sivity  ratio.  The  latter  ratio  is  dependent  on  the  electron  density 
because  the  two  lines  have  slightly  different  critical  densities. 
Using  our  value  for  the  electron  density  ~  100  cm-3  and  Fig.  3 
of  Lester  et  al.  (1987)  we  estimate  that  the  N/O  abundance  ratio 
to  be  ~  0.2  in  Cen  A.  This  value  is  consistent  with  the  range  of  ~ 
0.2  -  0.3  found  for  Galactic  HII  regions  (Rubin  et  al.  1988).  The 
nitrogen  to  oxygen  abundance  ratio  is  a  measure  of  the  chemical 
evolution  and  we  expect  it  to  increase  with  time  (cf.  the  solar 
value  of  ~  0.12). 


3.3.  Neutral  gas  lines 

Carbon  has  a  low  ionization  potential  (11.4  eV),  which  is  less 
than  that  of  hydrogen.  [Cil]  157  pm  line  emission  is  there¬ 
fore  observed  from  both  neutral  and  ionized  hydrogen  clouds. 
We  model  the  [C II]  line  emission  with  three  components:  Pho¬ 
todissociation  regions  (PDRs)  on  the  surfaces  of  UV  exposed 
molecular  clouds;  cold  (T  ~  100  K)  HI  clouds  (i.e.  the  cold 
neutral  medium  (CNM)  Kulkarni  &  Heiles  1987,  Wolfire  et  al. 
1995);  and  diffuse  HII  regions  (i.e.  the  warm  ionized  medium 
(WIM)  Heiles  1994). 


3.3.1.  HI  clouds 


It  can  be  shown  that  the  intensity  in  the  [C II]  line  emitted  from 
gas  clouds  with  density,  n(H)  and  temperature  (T)  is  given  by 
(c.f.  Madden  et  al.  1993) 


=  2.35  xlO"21  ' 


2  exp(=^) 
l  +  2exp{=^)  + 


—]XC+N(HI)  (1) 

nH 


where  the  critical  density  for  collisions  with  H,  ncrit  ~  3000 
cm'3  (Launay  &  Roueff  1977)  and  the  fractional  C+  relative 
to  hydrogen  is  Xc+  ~  Xc  =  1.4  x  10"4  (Sofia  et  al.  1997). 
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Fig.  3.  PDR  Lines  (Flux  Density  xlO17  W 
cm-2  /tm-1).  Note  the  [C  il]  lines  are  the 
total  observed  flux  density  i.e.  PDR  +  CNM 
+  WIM 


N(HI)  is  estimated  from  the  HI  21cm  map  of  van  Gorkom 
et  al.  (1990)  to  be  18.8  xlO20  atom  cm-2  at  the  SE  position. 
The  central  and  NW  positions  are  difficult  to  estimate  due  to 
HI  absorption  against  the  nuclear  continuum.  There  may  be  a 
central  hole  in  the  HI  and  the  column  density  is  certainly  not 
higher  than  the  peak  observed  in  the  SE  region  of  the  dust  lane 
(Van  Gorkom  et  al.  1990) 

Assuming  typical  Galactic  values  for  the  temperature,  T  ~ 
80  K,  and  hydrogen  density,  n  ~  30  cm'3,  results  in  an  estimated 
[C  II]  flux  of  3.5  x  10"20  W  cm"2  in  a  70  "  LWS  beam  at  the  SE 
position.  This  corresponds  to  4%,  1%  and  1%  of  the  observed 
[C  II]  flux  at  positions  SE,  Centre  and  NW  respectively.  The 
peak  HI  emission  line  flux  corresponds  to  1.9  x  10-19  W  cm-2 
which  is  only  6%  and  8%  of  the  [C  il]  flux  at  the  centre  and 
NW  positions  respectively.  We  conclude  that  there  is  very  little 
[C  II]  emission  in  our  beams  from  HI  clouds. 

3.3.2.  Diffuse  HII  regions 

Ionized  carbon  can  be  found  in  both  neutral  gas  and  ionized 
gas  clouds,  and  is  an  important  coolant  for  each.  We  detected 
[O  ill]  88  /tm  in  all  3  beam  positions  so  there  is  an  ionized 
gas  component  in  each  beam.  Using  the  constant  density  HII 
region  model  of  Rubin  (1985)  with  the  Kurucz  abundances, 
1049  ionizing  photons  per  second  and  our  derived  density,  ne  ~ 
100  cm-3  and  effective  temperature,  Teg  ~  35  500  K  we  can 
estimate  the  [C  II]  emission  from  the  HII  regions.  Applying  the 
model  [O  ill]  88  /rm  /  [C  II]  158  /im  line  ratio  of  0.35  to  the 
observed  [O  ill]  88  /im  line  flux  at  each  position  results  in  ~ 
10%  contribution  to  the  observed  [C  n]  line  flux  in  each  beam. 
Scaling  the  model  fluxes  to  the  distance  of  Cen  A  gives  ~  3000 
HII  regions  in  the  central  and  NW  regions  and  ~  1000  HII 
regions  in  the  SE  region. 


The  estimate  above  assumes  that  the  observed  lines  have  the 
same  filling  factor  in  the  large  LWS  beam.  If,  alternatively,  we 
were  to  assume  that  the  ionized  component  was  instead  dom¬ 
inated  by  a  contribution  from  an  extended  low  density  warm 
ionized  medium  (ELDWIM)  with  ne  ~  3  cm-3,  then  the  [C  il] 
flux  can  be  estimated  from  the  ratio  of  the  [C  n]/[N  il]  lines 
to  be  ~  18%  at  the  central  position.  The  observations  of  the 
[N  II]  121.9  /rm  line  at  the  NW  and  SE  positions  (with  lower 
signal  to  noise)  indicate  a  similar  fractional  component  (21% 
and  <  56%,  respectively). 

We  have  estimated  the  density  in  the  HII  regions  in  the  centre 
of  Cen  A  to  be  ~  100  cm'3  with  an  effective  temperature, 
T„ff  ~  35  500  K.  Based  on  the  HII  region  models  of  Rubin 
(1985)  we  estimate  that  ~  10%  of  the  observed  [C  II]  arises  in 
the  WIM. 

3.3.3.  PDRs 

Far-UV  photons  (6  eV  <  h//  <  13.6  eV)  from  either  O/B  stars 
or  an  AGN  will  photo-dissociate  H2  and  CO  molecules  and 
photo-ionize  elements  with  ionization  potentials  less  than  the 
Lyman  limit  (e.g.  C+  ionization  potential  =  1 1.26  eV).  The  gas 
heating  in  these  photodissociation  regions  (PDRs)  is  dominated 
by  electrons  ejected  from  grains  due  to  the  photoelectric  ef¬ 
fect.  Gas  cooling  is  dominated  by  the  emission  of  [0 1]  63  /tm 
and  [Cll]  158  /tm  emission.  Observations  of  these  lines,  the 
[O  1]  146  /tm  and  CO  (J=l-0)  2.6  mm  lines  and  the  FIR  con¬ 
tinuum  can  be  used  to  model  the  average  physical  properties  of 
the  neutral  interstellar  medium  (Wolfire  et  al.  1990).  Kaufman 
et  al.  (1999)  have  computed  PDR  models  over  a  wide  range  of 
physical  conditions.  The  new  code  accounts  for  gas  heating  by 
small  grains/PAHs  and  large  molecules,  and  uses  a  lower,  gas 
phase  carbon  abundance  (Xq  =  1 .4xl0-4,  Sofia  et  al.  1997)  and 
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Table  3.  PDR  diagnostic  line  intensity  ratios 


Line  Ratio 

Centre 

NW 

SE 

[O  i]  63  /xm  /  [C  n]  158  /xm 

0.7 

0.4 

0.6 

[O  i]  146  /xm  /  [O  i]  63  /xm 

0.06 

0.09 

<0.06 

([O  I]  63  /xm  +  [C  II]  158  /xm)/FIR 

0.006 

0.005 

0.006 

oxygen  abundance  (Xo  =  3.0xl0-4,  Meyer  et  al.  1998).  The 
[O  I]  63  ftm  /  [C II]  158  /it n  line  ratio  and  either  the  [O  i]  146 
/xm  /  [O  i]  63  gm  line  ratio  or  the  ([O  i]  63  /an  +  [C  II]  158  /itn) 
/  FIR  continuum  can  be  used  as  PDR  diagnostics  to  determine 
the  average  gas  density  (i.e  the  proton  density,  n  cm"  3),  the  av¬ 
erage  incident  far-UV  flux  (in  units  of  the  Milky  Way  flux,  Gc 
=  1.6  x  10-3  erg  cm-2  s-1)  and  the  gas  temperature. 

We  assume  that  the  measured  [C II]  flux  at  each  position 
should  have  ~  10%  subtracted,  due  to  the  HI  and  WIM  com¬ 
ponents,  before  it  is  used  to  model  the  PDRs  (if,  alternatively, 
a  20%  ELDWIM  contribution  is  subtracted  it  would  not  sig¬ 
nificantly  affect  the  PDR  parameters  derived  below).  The  PDR 
lines  are  plotted  in  Fig.  3  and  the  line  intensity  ratios  are  given 
in  Table  3. 

The  results  for  the  three  regions  are  consistent  with  each 
other,  having  a  gas  density,  n  ~  103  cm-3,  and  an  incident 
far-UV  field,  G  ~  102. 

At  the  NW  position,  only  the  combination  of  the  [O  I]  63 
/xm  /  [C  II]  158  /xm  ratio  and  the  ([O  I]  63  /xm  +  [C  II]  158  /xm) 
/FIR  continuum  ratio  gives  a  meaningful  solution  for  G  and  n. 
The  [O  I]  146  /xm  line  is  clearly  detected  but  with  a  very  rippled 
baseline  due  to  channel  fringes.  The  observed  [O  i]  146  /xm  line 
flux  would  need  to  be  reduced  by  ~  60%  in  order  to  obtain  a 
consistent  result  with  the  [O  I]  146  /tm  /  [O  I]  63  /xm  line  ratio 
predicted  by  the  PDR  model. 

The  LWS  results  for  the  nucleus  confirm  those  previously 
derived  from  IR,  submm  and  CO  observations.  The  consistent 
set  of  derived  PDR  conditions  for  all  three  positions  suggest 
that  the  observed  FIR  emission  in  a  70  "  beam  centred  on  the 
nucleus  is  dominated  by  star  formation  and  not  AGN  activity. 
Joy  et  al.  (1988)  mapped  Cen  A  at  50  and  100  /xm  on  the  KAO. 
They  concluded  that  the  extended  FIR  emission  was  from  dust 
grains  heated  by  massive  young  stars  distributed  throughout 
the  dust  lane,  not  the  compact  nucleus.  Hawarden  et  al.  (1993) 
mapped  Cen  A  at  800  /xm  and  450  /xm  with  a  resolution  of  ~10 
" .  They  attribute  the  large  scale  800  /xm  emission  to  thermal 
emission  from  regions  of  star  formation  embedded  in  the  dust 
lane.  They  note  that  the  H2  emission  within  a  few  arcseconds 
of  the  nucleus,  observed  by  Israel  et  al.  (1990),  indicates  that 
significant  UV  radiation  from  the  nucleus  does  not  reach  large 
radii  in  the  plane  of  the  dust  lane  i.e.  the  nuclear  contribution  to 
exciting  the  extended  gas  and  dust  disk  is  small. 

Eckart  et  al.  (1990)  and  Wild  et  al.  (1997)  mapped  Cen  A  in 
12CO  J=  1-0,  12CO  J=2-l  and  13CO  J=l-0.  All  three  maps  have 
two  peaks  separated  by  ~  90  "centred  on  the  nucleus.  It  is  inter¬ 
esting  to  note  that  our  SE  position  only  clips  the  lowest  contours 
of  the  CO  (1-0)  and  CO  (2-1)  maps  of  Wild  et  al.  (1997).  In  spite 
of  this  the  derived  PDR  parameters  are  consistent  with  those  en¬ 


compassing  the  bulk  of  the  molecular  emission.  There  must  be 
extended  low  level  CO  (1-0)  emission  beyond  the  sensitivity 
limits  of  the  Wild  et  al.  (1997)  maps.  The  lowest  contour  is  17.5 
Kkms-1,  corresponding  to  Mh2  ~  10s  M0  if  the  material  filled 
the  LWS  beam. 


4.  Summary  and  conclusions 

We  present  the  first  full  FIR  spectrum  from  43  -  196.7  /xm  of 
Cen  A.  We  detect  seven  fine  structure  lines  (see  Table  2),  the 
strongest  being  those  generated  in  PDRs.  At  the  central  position, 
the  total  flux  in  the  far-infrared  lines  is  ~  1%  of  the  total  FIR 
luminosity  (L43_ig7Mrra  =  3.2  x  109  L0  for  a  distance  of  3.5 
Mpc).  The  [C  II]  line  flux  is  ~0.4%  FIR  and  the  [O  I]  line  flux 
is  ~  0.2%  FIR.  These  are  typical  values  for  starburst  galaxies 
(Lord  et  al.  1996).  The  [O  ill]  52  /xm  /  [O  ill]  88  /xm  line  in¬ 
tensity  ratio  is  ~  0.9,  which  corresponds  to  an  electron  density, 
nc  ~  100  cm-3  (Rubin  et  al.  1994).  The  thermal  pressure  of 
the  ionized  medium  in  the  Cen  A  dust  lane  is  closer  to  that  of 
starburst  galaxies  (ne  ~  250  cm-3  in  M82  (Colbert  et  al.  1999) 
and  M83  (Stacey  et  al.  1999))  than  that  of  the  Milky  Way  (ne  ~ 
3  cm-3  (Petuchowski  &  Bennett  1993)). 

The  [N  ill]  /  [N II]  line  intensity  ratio  is  ~  1 .6,  giving  an 
abundance  ratio  N++/N+  ~  0.3,  which  corresponds  to  an  effec¬ 
tive  temperature,  Teff  ~  35  500  K  (Rubin  et  al.  1994).  Assuming 
a  coeval  starburst,  then  the  tip  of  the  main  sequence  is  headed  by 
08.5  stars,  and  the  starburst  is  ~  6  x  106  years  old.  If  the  burst  in 
Cen  A  was  triggered  by  the  spiral-elliptical  galaxy  merger  then 
its  occurance  was  very  recent.  Alternatively,  the  merger  trig¬ 
gered  a  series  of  bursts  of  star  formation  and  we  are  witnessing 
the  most  recent  activity. 

We  estimate  that  the  N/O  abundance  ratio  is  ~  0.2  in  the 
HII  regions  in  Cen  A.  This  value  is  consistent  with  the  range  of 
~  0.2  -  0.3  found  for  Galactic  HII  regions  (Rubin  et  al.  1988). 
N/O  is  a  measure  of  the  chemical  evolution  and  we  expect  it  to 
increase  with  time  (c.f.  the  solar  value  of  ~  0.12). 

We  estimate  that  ~  10%  of  the  observed  [C II]  arises  in  the 
WIM.  The  CNM  contributes  very  little  (<  5%)  [C  II]  emission 
at  our  beam  positions.  The  bulk  of  the  emission  is  from  the 
PDRs. 

We  derive  the  average  physical  conditions  for  the  PDRs  in 
Cen  A  for  the  first  time.  There  is  active  star  formation  throughout 
the  dust  lane  and  in  regions  beyond  the  bulk  of  the  molecular 
material.  The  FIR  emission  in  the  70  "  LWS  beam  at  the  nucleus 
is  dominated  by  emission  from  star  formation  rather  than  AGN 
activity.  On  scales  of  ~  1  kpc  the  average  physical  properties 
of  the  PDRs  are  modelled  with  a  gas  density,  n  --  103  cm-3,  an 
incident  far-UV  field,  G  ~  102  and  a  gas  temperature  of  ~  250 
K. 
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